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the Newtonian regime to a nonlocal and nonlinear modification of Poisson's equa- 
tion of Newtonian gravity. The nonlocally modified Poisson equation implies that 



CJ . nonlocality can simulate dark matter. Observational data regarding dark matter 

\^ \ provide limited information about the functional form of the reciprocal kernel, from 

which the original nonlocal kernel of the theory must be determined. We study this 
■ inverse problem of nonlocal gravity in the linear domain, where the applicability of 

. the Fourier transform method is critically examined and the conditions for the ex- 

Tlj" \ istence of the nonlocal kernel are discussed. This approach is illustrated via simple 

explicit examples for which the kernels are numerically evaluated. We then turn to a 
general discussion of the modified Poisson equation and present a formal solution of 
this equation via a successive approximation scheme. The treatment is specialized 
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\^ ■ to the gravitational potential of a point mass, where in the linear regime we recover 

the Tohline-Kuhn approach to modified gravity. 
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I. INTRODUCTION 

Newton's inverse-square force law for the universal attraction of gravity successfully ex- 
plained Kepler's laws of planetary motion and other solar-system observations. However, 
with the advent of Lorentz invariance, Einstein's general theory of relativity integrated 
Newtonian gravitation into a consistent relativistic framework that naturally explained the 
excess perihelion precession of Mercury and has thus far successfully withstood various ob- 
servational challenges. On the other hand, on galactic scales, the existence of dark matter 
has been assumed for decades in order to resolve observational problems, such as the "flat" 
rotation curves of spiral galaxies. As the nature of dark matter is still unknown, we take the 
view that what appears as dark matter may very well be a novel aspect of the gravitational 
interaction. Starting from first principles, a nonlocal generalization of Einstein's theory of 
gravitation has been developed in recent papers where it is demonstrated that non- 

local gravity can naturally simulate dark matter. However, the main motivation for the 
development of a nonlocal theory of gravitation came from the principle of equivalence, once 
the general outline of nonlocal special relativity had become clear [Sj. In nonlocal special 
relativity, nonlocality is induced by the acceleration of observers in Minkowski spacetime. 
Inertia and gravitation are intimately linked in accordance with the principle of equivalence 
of inertial and gravitational masses. Hence, gravity is expected to be nonlocal as well. That 
nonlocality can simulate dark matter emerged later in the course of studying the physical 



implications of nonlocal general relativity 



In this theory of nonlocal gravity, gravitation is described by a local field that satisfies 
nonlocal integro-differential equations. Thus gravity is in this sense fundamentally nonlocal 
and its nonlocality is introduced through a "constitutive" kernel. This nonlocal kernel is 
a basic ingredient of the gravitational interaction that must ultimately be determined from 
observation. A direct nonlocal generalization of the standard form of Einstein's general 



relativity would encounter severe difficulties js]; however, it turns out that one can also 
arrive at general relativity via a special case of the translational gauge theory of gravity, 
namely, the teleparallel equivalent of general relativity. It can be shown that this theory 
is amenable to generalization through the introduction of a constitutive kernel. In fact. 



the simplest case of a scalar constitutive kernel has been employed in [l-y] to develop a 
consistent nonlocal generalization of Einstein's theory of gravitation. In this approach to 
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nonlocal gravity, nonlocality survives even in the Newtonian regime and appears to provide 
a natural explanation for "dark matter" jl-lll; that is, nonlocal gravity simulates dark 



matter. Indeed, the nonlocal theory naturally contains the 



bhline-Kuhn phenomenological 



extension of Newtonian gravity to the galactic realm 12l-ll4 1 . 
Poisson's equation of Newtonian gravity, 

V^$jv(t,x) = 47rGp(t,x) 

is modified in the nonlocal theory to read js, 4| 

9k(x,y) 9<l>(y) 



(1) 



X 



E 



9x* 



d^y = A'KGp{y.) . 



(2) 



Here $jv is the Newtonian gravitational potential and any possible temporal dependence of 
the gravitational potential $ and matter density p has been suppressed in Eq. ([2]) for the 
sake of simplicity. Moreover, the nonlocal kernel k is a smooth function of u and f , so that 
k(x, y) = K{vL,v), where 

„:.^4M. (3) 



Thus V is the source of nonlinearity in Eq. ([2]). This nonlocal and nonlinear modification of 
Poisson's equation is invariant under a constant change of the potential, namely, $ t— ?■ $ + C, 
where C is a constant. In addition, Eq. (|2]) satisfies a scaling law; that is, if the matter 
density is scaled by a constant scale factor s, p ^ sp, then the potential is scaled by the 
same constant factor, $ i— t- s$. 

The functional form of the nonlocal kernel is unknown. Let us tentatively assume that it 
does have a dominant linear component A;(u) for some gravitational systems; that is. 



k = k{u.) + k(u, v) 



(4) 



where k{u, v) is a relatively small nonlinear perturbation. The physical justification for this 
supposition is that the implications of the linear form of Eq. ([2]), for the corresponding 
linear gravitational potential compare favorably with the physics of the fiat rotation 
curves of spiral galaxies . In the linear approximation, the scalar constitutive kernel k 
is only a function of u and we expect intuitively that even nonlocal gravity would weaken 
with increasing distance, so that k should go to zero as u := |u| — ?■ oo. Let us therefore 
start with Eq. (|2]) and consider, for simplicity, a linear kernel of the form k(x, y) := k{u), 
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so that 9k/9x* = —dk/dy^ in this case. Furthermore, let us assume that as y := |y| — )■ oo, 
|k(x, y)Vy$£(y) I falls off to zero faster than l/y^; then, using integration by parts and 
Gauss's theorem, Eq. ([2]) can be written in this case as 

V'$,(x) + y A;(x - y)V'$,(y)rf\ = 47rGp(x) . (5) 



We assume that this Fredholm integral equation of the second kind has a unique solution 15] , 
which can be expressed by means of the reciprocal convolution kernel g(u) as 

V'M^) = 47rGp(x) + AttG j g(x - y)p(y) A . (6) 

The conditions for the validity of this assumption will be examined in the following section. 
On the other hand, it is an immediate consequence of Eq. ([6]) that nonlocal gravity acts like 
dark matter; that is, the gravitational potential in the linear regime is due to the presence 
of matter of density p as well as "dark matter" of density p n given by 

pz5(x) = J g(x-y)p(y)rf=^y. (7) 

Thus the Laplacian of the gravitational potential is given by 47rG(p + p^), where, in the 
linear approximation, pd is the convolution of p and the reciprocal kernel q. In particular, 
for a point mass M, p(y) = M(5(y), we have po = Mq. 

A similar, but more intricate, result holds when nonlinearity is taken into consideration; 
that is, nonlocality still simulates dark matter, but the connection between p^ and p goes 
beyond Eq. ([7]). This assertion is based on the assumption that the nonlocal kernel is of the 
form of Eq. (jlj); then, Eq. ([2]) takes a similar form as Eq. ([5]), but with an extra source term 
iS(x) due to nonlinearity. That is, 

V2<l>(x) + y A;(x - y) V2$(y) A = 47rG[p(x) + 5(x)] . (8) 

Here iS(x) can be expressed as the divergence of a vector field, namely, 

5 = -V ■ i/, i/(x) = v)Vy^y)£y . (9) 

Assuming as before that Eq. (jH]) has a unique solution via the reciprocal convolution kernel 
g(u), we find 

V2$(x) = 47rG{p(x) + 5(x) + 1 g(x - y)[p(y) + S{y)]d'y} . (10) 
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It follows that the density of "dark matter" does have contributions from the nonlinear part 
of the nonlocal kernel. 

Neglecting nonlinearities, we note that the linear relationship between and p implies 
that one can write 

$,(x) = J g{^-y)p{y)d?y. (11) 
Here the Green function ^(u) is a solution of 

V^g{n) = A7rG[5{u) + g(u)] . (12) 

If p(x) = M5(x) for a point mass M in Eq. ([6]), we have $^ = MQ, which means that 
once the nonlocal gravitational potential is known for a point mass, one can determine the 
potential for any mass distribution via linearity as expressed by Eq. ffTTl) . Assuming that 
g(u) can be determined from observational data, the inverse problem must be solved to find 
the linear kernel k{u) from q{u). In this paper, we study this inverse problem in connection 
with the rotation curves of spiral galaxies; furthermore, we provide a general discussion of 
the solutions of Eq. ([2]) and investigate some of their physical implications. 

Consider the motion of stars within the disk of a spiral galaxy in accordance with the 
Kepler- Newton- Einstein tradition. For revolution on a circle of radius r := |x| about the 
central spherical bulge, the speed of rotation V of a star is given by V"^ = GAi/r, where Ai 
is the mass of the bulge, which we take to be the effective mass of the galaxy. Observational 
data indicate that is nearly constant; therefore, keeping the standard theory forces us 
to assume that mass Ai has a dark component that increases linearly with r. Assuming 
spherical symmetry for the distribution of this dark matter, we find that 

P«W = ii^. (13) 

where Vq is the constant asymptotic velocity of stars in the disk of the spiral galaxy. If 
dark matter does not exist, but what appears to be dark matter is in fact a manifestation of 
the nonlocal character of the gravitational interaction, then Eq. ([7]) together with Eq. f|T3|) 
implies that q{u) = M~^/(47rA) for p(x) = M5(x). Here M is the mass of the point source 
that represents the spiral galaxy and A = GM /Vq is a constant length. With this explicit 
form for the reciprocal kernel g, Eq. (E]) becomes identical to an equation that was first 
introduced by Kuhn in the phenomenological Tohline-Kuhn approach to modified gravity 
as an alternative to dark matter. 
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Let us digress here and briefly mention some relevant aspects of the linear phenomenol- 
ogy associated with the flat rotation curves of spiral galaxies. To avoid the necessity of 
introducing dark matter into astrophysics, Tohline 12| suggested in 1983 that the Newto- 
nian gravitational potential of a point mass M (representing, in effect, the mass contained 
in the nuclear bulge of a spiral galaxy) could instead be modified by a logarithmic term of 
the form 

^ , , GM GM ^ /r\ 

^^W = — + - . (14) 



r A VA. 

Here GM/X = Vq, where Vq is the approximately constant rotation velocity of stars in the 
disk of a spiral galaxy of mass M. Thus the constant length A is of the order of 1 kpc; hence- 
forth, we will assume for the sake of definiteness that A ~ 10 kpc. It follows that in this 
modified gravity model M oc Vq, which disagrees with the empirical Tully-Fisher law 16 |. 
The Tully-Fisher relation involves a correlation between the infrared luminosity of a spiral 
galaxy and the corresponding asymptotic rotation speed Vq. This relation, in combination 
with other observational data regarding mass-to-light ratio for spiral galaxies, roughly fa- 
vors M oc Vq, instead of M (x Vq that follows from Tohline's proposal. On the physical 
side, however, it should be clear that the Tully-Fisher empirical relation is based on the 
electromagnetic radiation emitted by galactic matter and thus contains much more physics 
than just the law of gravity for a point mass fisl. IitI. Tohline's suggestion was taken up and 
generalized several years later by Kuhn and his collaborators — see 13| and an illuminating 
review of the Tohline-Kuhn work in [l^. Indeed, in Kuhn's linear phenomenological scheme 
of modified gravity [l^, a nonlocal term is introduced into Poisson's equation, namely. 



1 f p(y) ,3. 



such that for a point source, p(x) = M5(x), Eq. (IT^ is a solution of Eq. ( TT^ . It follows 
immediately from a comparison of Eq. ([6]) with Eq. ( 1T5|) that 

Therefore, to make contact with observational data regarding the rotation curves of spiral 
galaxies, we suppose, as in previous work [l-4|, that the reciprocal kernel q in Eq. is 
approximately given by the Kuhn kernel in Eq. (fT6|) from the bulge out to the edge of spiral 
galaxies. 
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A remark is in order here regarding the nature of A. While for a point mass, A in Eq. ([T] 
is a universal constant, the situation may be different for the interior potential of a bounded 
distribution of matter. Consider, for instance, the Newtonian gravitational acceleration for 
a uniform spherical distribution of density po and radius R. The exterior acceleration has 
the universal form d^N/dr = GM/r"^ for r > R, where M = (47rpo/3)-R'^; as expected, this is 
identical to the acceleration for a point source of fixed mass M. For < r < R, the interior 
Newtonian gravitational acceleration at a fixed radius r is given by d^^/dr = [GM/R^)r, 
which decreases with increasing R when M and r are held fixed. Extrapolating from this 
natural consequence of Newtonian gravity to the nonlocal domain, one might expect that 
in the interior of spiral galaxies, for instance, A in Eq. ( fT5l) might depend on the size of 
the system. Indeed, in Kuhn's work. A, with a magnitude of more or less around 10 kpc in 



HQ. 



Eq. f[T5|) . is taken to be larger for larger systems 

Rotation curves of spiral galaxies can thus provide some information regarding the nature 
of the reciprocal kernel g(u) in the linear case, where Eq. ([6]) can be directly compared with 
the Tohline-Kuhn scheme. However, to determine the corresponding nonlocal kernel k{u), 
we need to know the functional form of g(u) over all space. The simplest possibility would 
be to assume that q(u) = u~^/(47rA) holds over all space; however, the corresponding k{u) 
does not exist — see [2|, especially Appendix E, for a detailed discussion of this point. We 
therefore take up the crucial question of the existence of the linear kernel k{u) for spiral 
galaxies in sections II and III. The general nonlinear inverse problem of nonlocal gravity is 
beyond the scope of this paper; therefore, we concentrate in sections II and III on the linear 
problem of finding k{u) from certain very simple extensions of q{u) beyond q{u) oc that 
is implied (in the galactic disk) by the flat rotation curves of spiral galaxies. We then turn 
to the other main goal of this paper, which is to tackle the general nonlinear form of Eq. ([2]). 
Unfortunately, the general form of the nonlocal kernel is unknown at present. Nevertheless, 
a formal treatment of Eq. ([2]) is developed in Sec. IV without recourse to Eq. (jlj) or any 
specific assumption about the nature of the nonlocal kernel. As an application of this new 
approach, we consider the gravitational potential of a point mass in Sec. V. Physically, we 
regard the point mass as an idealization for a spherical distribution of matter; in practice, it 
will stand, for instance, for the mass of a spiral galaxy, most of which is usually concentrated 
in a central spherical bulge — see, in this connection, section IV of j^. The linear regime, 
where the kernel is independent of v, will be investigated in the first subsection of Sec. V; in 
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fact, we illustrate the effectiveness of our procedure by demonstrating how previous results 
can be recovered in the new setting. Specifically, we show that the Tohline-Kuhn scheme 
can be recovered in this case from the general treatment of Sec. IV. Finally, Sec. VI contains 
a discussion of our results. 

II. INVERSE PROBLEMS IN LINEAR NONLOCAL GRAVITY 

The nonlocal theory of gravitation under consideration in this paper is based on the 
existence of a suitable "constitutive" kernel k. For the gravitational potential of spiral 
galaxies, we assume at the outset that k contains a dominant linear part k{u) and a nonlinear 
part that we can ignore in the context of the present discussion. The reciprocal of k{u) for a 
spiral galaxy is expected to be of the form of the Kuhn kernel f|T6|l in order that the nonlocal 
theory could simulate dark matter and be therefore consistent with observational data. It 
turns out that with the simple form of q given in Eq. ( |T6l) . the corresponding reciprocal 
kernel does not exist. That is, if we naturally extend the simple Kuhn kernel to the whole 
space beyond a galaxy, then there is an infinite amount of simulated dark matter and we 
have a q{u) oc for which there is no finite k{u). However, the nonlocal theory is based 
on the existence of a finite smooth nonlocal kernel. This important problem is taken up in 
this section. To determine k{u) from its reciprocal, it is necessary to extend the functional 
form of kernel ( IT6l) so that it becomes smoothly applicable over all space and falls off rapidly 
to zero at infinity. Let q{u) be this extended reciprocal kernel. We must ensure that its 
reciprocal k{u), the constitutive kernel of nonlocal gravity in the linear regime, indeed exists 
and properly falls off to zero as m — !■ oo. To this end, we show in this section that it is 
sufficient to require that q{u) and k{u) be smooth absolutely integrable as well as square 
integrable functions over all space, so that their Fourier integral transforms exist as well. 

In the linear Newtonian regime of nonlocal gravity, the nonlocally modified Poisson 
equation is a Fredholm integral equation of the second kind, cf. Eq. (j5]). 



which we assume has a unique solution that is expressible by means of the reciprocal kernel 

q as 




(17) 




(18) 
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The inverse problem, however, involves finding k{u) once g(u) is known; that is, we wish to 
obtain Eq. f lT7|) starting from Eq. f|T8|) . Let us note that our methods can be used in either 
direction due to the obvious symmetry of Eqs. ( |T71) and (ITSl) . 

It is useful to express Eq. ( ITTl) in operator form as (/ + IC)g = /, where / is the identity 
operator and /C is the convolution operator JCg = k-k g. Formally, we expect the solution of 
this equation is (/ + Q)f = g, where Qf = q f- Moreover, f = {I + ]C)g = (/ + Q)~^g 
would be equivalent to Eq. f lT7|) . From a comparison of Eq. f|T7|) with Eq. ([5]), we see that the 
quantities of interest are / = AnGp and g = V^$f. Here, the function p models the density 
of matter in space and is the linear gravitational potential. Both of these functions can be 
considered to be smooth in the continuum limit for matter distributions under consideration 
throughout this work. 



A. Liouville-Neumann Method 



It is formally possible to obtain Eq. (fTTll from Eq. ( ITSll . or the other way around, by the 



application of the Liouville-Neumann method of successive substitutions 15|. That is, we 
start with Eq. f|T8|) and write it in the form 

/(x) = (7(x)- J g(x-y)/(y)rf\. (19) 

Then, we replace /(y) in the integrand by its value given by Eq. f lT9|) . Iterating this process 
eventually results in an infinite series — namely, the Neumann series — that may or may not 
converge. 

A uniformly convergent Neumann series leads to a unique solution of the Fredholm in- 
tegral equation ( !T9|) . Moreover, one can determine k in terms of q. The procedure for 
calculating kernel k of Eq. f lT7|) in terms of the iterated q kernels has been discussed, for 



instance, in jsj; however, a sign error there in the formula for iterated kernels must be 
corrected: An overall minus sign is missing on the right-hand side of Eq. (3) of jsl. The 
spherical symmetry of q, in the cases under consideration in this section, implies that all 
iterated kernels are functions of u. Thus let qniu), n = 1,2,..., be the relevant iterated 
kernels such that qi = q and 

g„+i(|x -y\) = -J g(|x - z|)g„(|z - y\)dh . (20) 
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It follows that in this case the nonlocal kernel is 

oo 

M|x-y|) = -5^gn(|x-y|). (21) 

n=l 

In trying to determine k from q, one can therefore start from the study of the Neumann 



series. Using the approach developed in 15|, and working in the Hilbert space of square- 



integrable functions, it can be shown by means of the Schwarz inequality that the Neumann 
series converges and the nonlocal kernel exists for ||g|| < 1; moreover, the solution of the 
Fredholm integral equation ( ITSl) by means of the Neumann series is unique. However, the 
norms of the convolution operators of interest in our work are not square integrable, since 

||g||2 = j q'^{x-y)£xd^y = j d^x j q'^{u)d^u = oo , (22) 



so that the standard Hilbert space theory developed in 15[ is not applicable here. That is, 
the norm of q could be finite, but the norm of the corresponding convolution operator in 
is infinite. 

On the other hand, let us suppose that the space of functions of interest is a Banach 
space B and that Q is a bounded linear operator on i3, Q : B ^ B, with ||Q|| < 1; then, 
one can show that (/ + Q) has an inverse given by Xl^o(~2)"' "where the series converges 
uniformly in the set of all bounded linear operators on B. In this formula, (— Q)" in the 
series corresponds to the iterated kernel g„, n = 1,2,3, in the Neumann series and the 
convergence of the series is equivalent to the existence of kernel k. It seems that the sufficient 
condition ||Q|| < 1 for the convergence of the Neumann series cannot be satisfied for any 
physically reasonable extension of the Kuhn kernel ( !T6|) : that is, our various attempts in 
this direction have been unsuccessful. In short, the Neumann series diverges; therefore, we 
resort to the Fourier transform method in this work. 

We caution that our mathematical approach may not be unique, as the theory may 
work in other function spaces that we have not considered here. The general mathematical 
problem is, of course, beyond the scope of this paper. 



B. Fourier Transform Method 



Following well-known mathematics (see, for example, 18|, §9.6]), precise conditions can be 
determined for our convolution operators to be invertible. The basic idea is to use the Fourier 
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transform T defined for functions (that is, functions which are absolutely integrable over 
all of space) by 



to prove a basic lemma: Ij q is in , then its Fourier transform q is continuous and the 
convolution operator Q given by Qf = q-kf is a bounded operator on whose spectrum is the 
closure of the range of the q. The proof outlined here uses standard results of theory. By 
an elementary argument, the Fourier transform of an arbitrary function is continuous. A 
deeper result (Plancherel's theorem) states that the Fourier transform can be extended to a 
bounded invertible operator on that preserves the norm; in other words, the extended 
Fourier transformation is an isometric isomorphism of the Hilbert space L^. This extended 
operator (also denoted by J-" and simply called the Fourier tranform) maps convolutions to 
products: J^[f g] = J-'[f]J-'[g]. Let Q be the (multiplication) operator defined on by 
Qf = qf, where q is the Fourier transform of the kernel of Q. If follows from the definitions 
and the action of the Fourier transform on convolutions that Q = T^^Q,T\ therefore, the 
spectra of the operators Q, and Q, coincide. The spectrum of the multiplication operator Q 
(with its continuous multiplier q) on is the closure of the range of g, as required. 

A corollary of the lemma is the result that we require to analyze our integral equations: 
// the number —1 is not in the closure of the range of q, then I + Q is invertible. 

To see how these ideas can be used to obtain Eq. f|T7j) more explicitly, we consider 
integral equation (fT8l) in the form f + q-^ f = g and apply the Fourier transform — under the 
assumption that / and g are in and q is in L} — to obtain the equivalent equation 




(23) 



{l + q)f = 9- 



(24) 



If 1 + g ^ 0, then 




-)g = g + kg, 



(25) 




(26) 



Applying the inverse Fourier transform, it follows that 



f = g + J^-'[kg]. 



(27) 



// there is an L 



function k such that J^[k] = k, then 



f = 9 + k-kg 



(28) 
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and I + JC, where JCg = k-k g, is the inverse of I + Q. 

These results are illustrated in detail in the next section via specific examples. In these 
applications, we will employ a useful lemma regarding the Fourier sine transform. Consider 
the integral 

POO 

J(0 = / h{x,Osm{^x)dx. (29) 
Jo 

For each C, € (0,oo), let h{x,$,) be a smooth positive integrable function that monotonically 
decreases over the interval of integration; then, J(^) > 0. To prove this assertion for each 
> 0, we divide the integration interval in Eq. fl29|) into segments {2'K^~'^n, 2Ti^'~^n + 2Tx^~^) 
for n = 0,1,2,.... In each such segment, the corresponding sine function, sin(,^x), goes 
through a complete cycle and is positive in the first half and negative in the second half. 
On the other hand, the monotonically decreasing function h{x,^) > is consistently larger 
in the first half of the cycle than in the second half; therefore, the result of the integration 
over each full cycle is positive and consequently J{^) > 0. 

For ^ — )■ 0, sin(^x) — and hence J(0) = 0, while for ^ — j- oo, the integration segments 
shrink to zero and J tends to in the limit as ^ — )■ oo, if the corresponding limit of h{x,^) 
is finite everywhere over the integration domain. This latter conclusion is, of course, a 
variation on the Riemann-Lebesgue lemma. 



III. EXISTENCE OF THE LINEAR KERNEL: EXAMPLES 

We wish to determine k{u) from a knowledge of q{u). Let us note that if g(u) is only a 
function of the radial coordinate u, Eq. ( p3|) reduces to 

g(^) = — / rq{r) sm{^r)dr , (30) 
s Jo 

where ^ : — oo is the magnitude of ^. That is, in Eq. (123|) . we introduce spherical polar 
coordinates (r, 6, (p) and imagine that the coordinate system is so oriented that $, points 
along the polar axis; then, the angular integrations can be simply carried out using the fact 
that 



In general, the Fourier transform g(^) of a square-integrable function g(u) is square inte- 
grable. In case 1 + g 7^ and k{^) is an function, the Fourier transform method of the 
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previous section becomes applicable here. Thus a suitable nonlocal kernel k{x.) exists in this 
case and is given by 



For the explicit radial examples under consideration, Eq. fl32|) takes the form 

A. An Example 

To extend Kuhn's kernel, Eq. f lT6|) . smoothly over all space, one may consider, for instance. 



1 d 



47rA du 



F(u)i 



a + u. 



^(«) = ^3 — ■ (34) 



Here a is a constant length scale characteristic of the radius of nuclei of spiral galaxies, 
while F{u) is a smooth function that is nearly unity over much of the interval (0,^4), but 
then rapidly drops off to zero as m — t- oo. The constant A represents another length scale 
characteristic of the radius of galactic disks. Thus a X < A under physically reasonable 
conditions, and for a u < A, one can show that Eq. (!34|) essentially coincides with 
Kuhn's kernel. To recover the flat rotation curves of spiral galaxies, Eq. f lM|) should agree 
with Kuhn's kernel from the bulge, which is, say, at a distance of about A ~ 10 kpc from the 
galactic center, out to the edge of the galaxy, which is, say, at a distance of about 3A. The 
function F{u) can be chosen so as to render 1 + g > 0; to illustrate this point, we choose 

F{u) = e-"/^ . (35) 



Let us now consider the application of Eqs. (1301) and (1331) to the particular case of Eqs. 
- (l35l) . Regarding the parameters that appear in these equations, it is sufficient to suppose 
at the outset only that A > 0, a > and A > are length scales of interest; moreover, we 
introduce, for the sake of simplicity, 

a:=j. (36) 

We note that q{u) is smooth and positive everywhere and rapidly decreases to zero at 
infinity; indeed, q{u) is integrable as well as square integrable. It is preferable to work with 
dimensionless quantities; thus, we let all lengths — such as r, u, A, a and A — be expressed in 
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units of A. Then, X^q, A^, Aa and q are dimensionless. Similarly, X^k and k are dimensionless. 
Henceforth, we deal with these dimensionless quantities; in effect, this means that A = 1 in 
the following formulas. It is possible to show that in this case for any C, > 0, q{C,) > —a. 
Substituting Eqs. (134|) and (135|) in Eq. ( 130|) and integrating by parts, we find 

m = -J^ ^-[rsin(er)]cir. (37) 

Next, differentiating rsin(^r) and noting that 

sin(,^r) + cos(.^r) = [sin(,^r) — cos(,^r)] + (a + r)^ cos(,^r) , (38) 

Eq. (I57|) can be written as 

POO 

g(^)=X+/ e""" COS (^r)dr, (39) 
Jo 



where 



r e~^^ cos{^r)dr = , (40) 

Jo ,— , " + r 



according to formulas 3.893 on page 477 of Ref. 

roo ^—ar 



19|, and 



sin(^r) — cos(^r)] dr . (41) 



Uo « + ^ 

Let us now introduce an angle 7 connected with such that 

:= tan7 , (42) 

and note that as ^ : — t- 00, we have 7:0—7- 7r/2 and 7/^ : a — t- 0. It is useful to introduce 
a new variable X in Eq. (HT|) . r = X + 7/^, since 

sin(fr) — af cos(fr) = — ^ — sinffr — 7) . (43) 
^ C0S7 

Then, Eq. ( 14T]) can be written as 

I = / -—- sin(eX) dX . (44) 

ecos7y_^/^ a + l + X ' ^ ' 

In this expression, the integration from X = —7/.^ to 00 can be expressed as a sum of two 
terms, one from to 00 and the other from X = —7/^ to 0. That is, X = P + A/", where 

-07/5 poo 

V = / sin(eX) dX (45) 

4 cos 7 Jo a + j + X 
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is positive by the argument presented at the end of Sec. II, since exp(— aX)/(X + + 7/^) 
is a smooth positive integrable function that monotonically decreases for X : — 00, and 

-07/5 rO -aX 

M = / sin(eX) dX , (46) 

Ccos7 J_^/^ a + ^ + X 

which is negative. This latter point can be made explicit by introducing a new variable 
= —Y into Eq. ( l46l) : then, we have 

- e cos 7 AfiO = / 7 TTV— T ^ • (-^7) 

The right-hand side of this equation involves an integrand that increases monotonically from 
to sin 7/(0^) as y : — 7. Thus the right-hand side of Eq. (147|) is less than 7 sin 7/(0^); 
consequently, A/'(0 > -7/^ by Eq. (j42]). As < 7/^ < a, we conclude that A^(0 > -a. 
Collecting our results, we therefore have that q > —a and 

l + g(0>l-a. (48) 



Hence the Fourier transform method of the previous section is applicable to Eq. ( 1341) if a < 1. 
It then follows from Eq. ( I32l) that \k\ < |g|/(l — a), so that k{^) is in as well, and we can 
find the nonlocal kernel k from Eq. fl33|) . 

In connection with the rotation curves of spiral galaxies, it is useful to consider the 
amount of sham dark matter that is associated with such a model. An estimate of the net 
amount of simulated dark matter Md can be obtained from the integration of Eq. (I7l) over 
all space, where we set p(y) = M6{y) for the sake of simplicity. Thus 

f« 47r / r^q{r)dr . (49) 

^ 

For the example under consideration, we find from Eqs. ( l34l) . (135|) and ( l49l) via integration 
by parts that 

r (50) 



M a Jq aa + ( 

where the definite integral can be expressed in terms of the exponential integral function as 



1 + aa exp(aa)Ei(— aa) — see, for instance, page 311 of Ref. [19|. Thus nonlocality in this 
case simulates, in effect, a net amount of dark matter that is nearly 2/a times the actual 
mass M of the galaxy, since the integral term on the right-hand side of Eq. ( l50i) is nearly 
unity for physically reasonable values of the parameters, namely, < aa ^ 1. Indeed, for 
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< X ^ 1, Ei(—x) y C + \nx, where C = 0.577... is the Euler-Mascheroni constant — see 



19|. Further considerations involving kernels q and k are relegated to 



page 927 of Ref. 
Appendix A. 

Choosing dimensionless parameters a = 0.1 and a = 0.001 in this example, the corre- 
sponding numerical results for q and k are presented in Figures 1 and 2. 

B. A Second Example 

The purpose of this subsection is to discuss a second example of a nonlocal kernel. We 
start with the reciprocal kernel 

,(„)^-L l + °(°+;) e-". (51) 

Here a := 1/A, as before. This reciprocal kernel has a central cusp and behaves as 1/u for 
u <^ a, which is reminiscent of the density of dark matter in certain dark matter models — 



see, for instance, Ref. {q]. Moreover, for a <^ u < A, q behaves like the Kuhn kernel, while 
for u ^ A, it falls off exponentially to zero. 

Kernel ( 15T|) is a smooth positive integrable function that is in L^. Using dimensionless 
quantities, Eq. fl30|) takes the form 

1 f°° 1 

q(n = - (a + )e-"'^ sin(^r)dr . (52) 

' J1 + r 



From formulas 3.893 on page 477 of Ref. 



19| . we find 



-1 POO —ar 
a i / 6 

««) = ^ + J„ — smKO*. (53) 

Here we can directly use the lemma given at the end of Sec. II, since exp(— ar)/(a + r) is a 
smooth positive integrable function that decreases monotonically for r : — oo, to conclude 
that for < < oo, q[^) > 0, while q{C,) — )■ as — oo by the Riemann-Lebesgue lemma. 
It follows that \k{C,)\ < 1^(01' so that k is in as well and the nonlocal kernel k can be 
determined via Eq. fl5^ . 

In this case, the analog of Eq. fH9|) is given by 

^ ^ - + ae""Ei(-aa) , (54) 

which is, for < aa ^ 1, nearly the same as in Eq. ( I50l) . For this second example, the 
numerical results involving q and k for a = 0.1 and a = 0.001 are presented in Figures 3 
and 4. 



17 



20 
15 
10 

5 



2 4 6 8 10 

FIG. 1: Plot of q versus A.^ for the reciprocal kernel q given in Eqs. ()34p and (1350 . The parameter 
values are Aa = 0.1 and a/A = 0.001. The function q starts from q[Qi) ~ 20 and rapidly falls off 
initially, but then slowly decreases to zero as A^ — > oo. 
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FIG. 2: Plot of — A^fc versus u/\ when the reciprocal kernel q is given by Eqs. (j34p and (j35p . The 
parameter values are Aa = 0.1 and a/A = 0.001, just as in Figure 1. The function — A'^A; starts 
from ~ 80101 at u = and drops off to nearly zero very fast; in fact, for n/A > 2.5 it is essentially 
zero at the level of accuracy of this plot. 

Appendix A contains further useful mathematical results relevant to the examples de- 
scribed in this section and the corresponding numerical work presented in the figures. 

The nonlocal "constitutive" kernel k{u) turns out to be negative for models of spiral 
galaxies under consideration in this work. Moreover, as Figures 2 and 4 indicate, nonlocality 
in this case involves sampling sufficiently close spatial regions. Indeed, around any point x, 
the influence of the field amplitude at point y may be significant only when u = |x — y| is 
smaller than around 2.5A, or about 25 kpc. In fact, as expected, at any given field point, the 
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FIG. 3: Plot of q versus for the reciprocal kernel q given in Eq. ()5ip . The parameter values 
are Xa = 0.1 and a/A = 0.001, just as in Figure 1. As pointed out in Appendix A, q is in this 
case always slightly larger than the one in Figure 1, but this is hardly noticeable for the parameter 
values under consideration. For instance, at A^ = 2, g w 0.773 in Figure 1, while here q ^ 0.779. 
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FIG. 4: Plot of —X^k versus u/X when the reciprocal kernel q is given by Eq. (j51|) . The parameter 
values are Xa = 0.1 and a/A = 0.001, just as in Figure 3. The function — A^A; starts from oo at 
u = and drops off to nearly zero very fast; in fact, for u/X > 2.5 it is essentially zero at the level 
of accuracy of this plot. Though this figure appears to be indistinguishable from Figure 2 in the 
plotted range, their numerical values are indeed different. 

nonlocal influence of the field amplitude at a nearby point decreases with increasing distance 
extremely fast. 
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IV. NONLOCAL AND NONLINEAR POISSON'S EQUATION 

The purpose of this section is to present the main outhnes of a formal approach to the 
modified Poisson equation with a general nonlinear kernel. That is, we do not assume here 
that the nonlocal kernel consists of a dominant linear part together with a small nonlinear 
perturbation. 

The right-hand side of Eq. ([2]) can be replaced by the Laplacian of the Newtonian potential 
via Eq. ([1]). Furthermore, it is straightforward to see that the nonlocal contribution to Eq. ([2]) 
can be written as the divergence of a vector field. It follows from these remarks that modified 
Poisson's equation can thus be written as V ■ = 0, where 

* = V$ + j k(x, y) Vy$(y)rf^y - V$7v • (55) 

For a bounded matter distribution, we can write the solution of Eq. ([T]) as 

$^(x) = -G f T^^4'y. (56) 

J y I 

so that, as is customary, the Newtonian gravitational potential is assumed to be zero at 



infinity 
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We are interested in the solution of the nonlinear integral equation 

V$ + y" k(x, y) Vy$(y)rf^y = + * . (57) 

The divergence-free vector field ^ must be such that Eq. ( 1571) is integrable. Indeed, the 
integrability condition for Eq. fl57j) is that 

j Vxk(x, y) X Vy^{y)(fy = V x * . (58) 

In effect, the modified Poisson equation has thus been once integrated and reduced to 
Eqs. fl57|l and f l58|l . which are, however, still unwieldy. 

Let U represent the left-hand side of Eq. (158|) : then, one can express Eq. (l58l) as V x ^' = 
U, where U is divergence- free. The divergence of ^ vanishes and its curl is U; therefore, 
the curl of — U is the Laplacian of 

V2* = -VxU. (59) 

If U(x) is bounded for small r = |x|, falls off to zero faster than 1/r for large r and ^' — ?■ 
as r — )■ oo, then 

*(x) = i-V X y ^Zjfy ' U(x) = V X y k(x, y) Vy$(y)c/3y . (60) 
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In this way, one can determine ^ from the integrabihty condition, namely, Eq. (158|) . and 
substitute it back into Eq. ([5j 



A. Formal Solution via Successive Approximations 

The solution of the integral equation for the gravitational potential $ is expected to 
consist of the Newtonian gravitational potential $ together with nonlocal corrections as in 
a Neumann series; therefore, it is natural to devise a formal solution of Eq. ( l57l) using the 



method of successive approximations 15[. In view of the divergence of the Neumann series in 
connection with the considerations of the first part of this paper, we must assume here that 
the gravitational system under consideration in this subsection cannot be approximated by a 
point mass. Let $07 ^n, ■■■ be a series of approximations to the gravitational potential 

such that $0 •= "^AT and $„ approaches $ in the limit as n — )■ 00. Moreover, for n > 0, we 
define to be such that 

V<l>i = V<^N + *o - y K{n, t;o)Vy$o(y)rf\, (61) 



= V^N + *n-i - j K{u, Vn-i)Vy^n-i{y)d'y, (62) 

V$„+i = V^N + ^n- j K{n, Vn)Vy<!>n{y)d'y, (63) 

and so on. Let us recall here that k;(x, y) = A'(u, v) and we have extended the definition of v 
in Eq. (|3]) such that Vn = |Vy$„(y)|/| Vx$n(x)|. Moreover, ^n, ■■■ are such that 

the integrabihty condition is satisfied at each step of the approximation process, namely, 

V X = J VxA'(u,w„) X Vy$„(y)rf3y (64) 

for n = 0,l,2,.... Here for instance, can be expressed in terms of $„ using the method 
described in Eqs. f l59|) and fl60|) : then, the result may be employed in the expression for 
V$„+i in Eq. (1^ of the successive approximation scheme. As n — )■ oo, we expect that 
approaches $ and approaches so that the limiting form of Eq. f lM|) coincides with 
Eq. ( I58l) . The convergence of this successive approximation process depends of course upon 
the nature of the kernel and its treatment is beyond the scope of this work, as the general 
form of kernel K{u,v) is unknown at present. 
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The general solution of Eq. (|2]) presented in this section can be used, in principle, to 
restrict the form of kernel K{u, v) on the basis of observational data. Nonlocal gravity 
simulates dark matter jl-^; therefore, it may be possible to determine the general nonlinear 
kernel K from the comparison of our general solution of Eq. ([2]) with astrophysical data 
regarding dark matter. However, the treatment of this general inverse problem of nonlocal 
gravity is a task for the future. 



V. GRAVITATIONAL POTENTIAL OF A POINT MASS 

As a simple application of the formal procedure developed in the previous section, we will 
consider here the gravitational potential due to a point mass M at the origin of spatial co- 
ordinates, so that p(x) = M(5(x). The corresponding Newtonian potential is $Ar = —GM/r 
and we expect that $ is also just a function of r as a consequence of the spherical symmetry 
of the point source. Similarly, it is natural to assume that the kernel's dependence on u is 
only through its magnitude u due to the isotropy of the source. A detailed investigation 
reveals that ^' = in this case; we outline below the main steps in this analysis. 

In computing the integral term in Eq. fl57|) . we introduce the spherical polar coordinate 
system {y, ^, if) in which the polar axis is taken to be along the x direction. The kernel in 
Eq. flFr|) is then just a function of r, y and cos^9; moreover, Vy$(y) equals d^{y)/dy times 
the unit vector in the y direction. The azimuthal components of this unit vector vanish upon 
integration over all angles and only its polar component remains. Therefore, ^ is purely 
radial in this case, namely, 

* = x(r)x, (65) 

which satisfies the integrability condition given in Eq. (158|1 . since in this case the curl of ^ 
identically vanishes. Here x(r) can be determined from the requirement that V ■ ^' = 0. It 
then follows that x = Tn/r^, where m is an integration constant. Thus ^ = m:K/r^; that is, 
the right-hand side of Eq. (!57|) is radial in direction and is given by {GM + m)x/r3. The 
resulting ^ in effect indicates the presence of an extra delta-function source at the origin of 
spatial coordinates. We therefore set m, which is effectively a new mass parameter, equal 
to zero, as it simply renormalizes the mass of the source. Thus ^' = and Eq. fl57|) reduces 
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in this case to 

^ + 27r r r K{u, v) ^^y^dy cos ^ sin M^ = ^, (66) 
dr Jq Jq dy r 

where u and v are given by 



u = Jr^ + y^ -2rycos^, = . (67) 

dy dr 

The extra factor of cosi? in Eq. (!66|) is due to the fact that in Eq. (l57|l the component of 
y/y along the polar axis is cos-i?. 



A. Linear Kernel 

Let us assume, for the sake of simplicity, that k(x, y) = k{u), so that in this subsection 
we are only concerned with a nonlocally modified Poisson's equation that is linear with a 
kernel that depends only on m as a result of the spherical symmetry of the point source. 
Then, Eq. fl^^ for the linear gravitational potential reduces to 

^ + 2tt [ I k(Jr^ + 2/2 _ 2ry cos ^ ) ^^$^y^dy cos d sin = ^ , (68) 
dr Jo Jo dy r^ 

which means that a linear integral operator with kernel k acting on d^g/dr results in GM/r"^. 
We note in passing that the successive approximation method of the previous section leads 
in this case to the standard Liouville-Neumann solution of Eq. ( 168|) via iterated kernels of 
the Fredholm integral equation of the second kind 15|; however, as discussed Sec. II, the 
Neumann series diverges in this case and the corresponding solution does not exist under 
physically reasonable conditions. Therefore, we adopt the Fourier transform method and 
let q{u) be the kernel that is reciprocal to k{u)] then, d^g/dr is given by the linear integral 
operator, with k replaced by q, acting on GM/r^. That is, 

d^e GM 



dr r 



2 



PTT POO 

+ 2ttGM / / q{ v^r^ + _ 2ry cos )dy cos d sin ddd . (69) 
Jo Jo 



Substituting the Kuhn kernel (fT6|) for q in Eq. (!69|) and performing the ^/-integration in 
the resulting integral first, we find that for d G (0, vr], 

dy 71 — 13 



(y — r cos + r"^ sin^ ■& rsin-t? 



(70) 



23 



The t^-integration is then straightforward and the end result is 

d^i GM GMl 
ar A r 

in agreement with the radial derivative of Eq. f[T^ . In this way, starting from our general 
solution of the modified Poisson equation, we again recover the Tohline-Kuhn scheme of 
modified gravity. 

B. Nonlinear Kernel 

To gain some insight into the role of nonlinearity in Eq. fl66p . let us suppose that nonlin- 
earity constitutes a very small perturbation on a background linear kernel. In fact, we set 
K{u,v) = k{u) + eP{u,vg), where e, < e ^ 1, is a sufficiently small parameter and Vi is 
obtained from v in Eq. f l67|) by replacing $ with We thus expand $ to first order in e 
and thereby develop a simple linear perturbation theory for Eq. (166|) such that 

$ = $^ + e^ne . (72) 

Here $^(r) is given in general by Eq. (169|) and is the perturbation potential due to 
nonlinearity. Moreover, Eq. f l66|) implies that 

^^ + 271 [ [ k{^r^ + y2 _ 2rycosd) '^^'f'^\ ^dycosdsiYidM = N{r) , (73) 
dr Jq Jq dy 

where N{r) is due to the nonlinear part of the kernel and is given by 

N{r) = -2tx f [ P{u, ve) ^"^/^^ y^dy cos ^ sin . (74) 
Jo Jo dy 

As in the previous subsection, Eq. ( TfSl) can be solved by means of kernel q{u) that is 
reciprocal to k{u) and we find 



TT POO 



^ = N{r) + 2n I I q{\/r'^ + — 2ry cos 1!} ) N {y)y'^ dy cos i9 sin Mi^ . (75) 
dr Jo Jo 

A consequence of this result should be noted here: Inspection of Eqs. fl69|) . (17^ and f l75]) 
reveals that ^neir) is simply proportional to the gravitational constant G. This feature is 
an example of the general scaling property of Eq. ([2]), which implies that any solution $ of 
Eq. ([2]) must be proportional to G. 

It is therefore possible to see that our nonlocal as well as nonlinear modification of New- 
tonian gravity cannot behave as in the Modified Newtonian Dynamics (MOND) approach 
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to the breakdown of Newtonian gravity 2ll423l|. From the scahng property of our modified 



Poisson's equation, we expect that the gravitational potential $ is in general proportional to 
the gravitational constant G, since the source term in Eq. ([2]) is proportional to G. There- 
fore, the nonlocal theory, as a consequence of its particular nonlinear form in the Newtonian 
domain, does not contain a MOND regime, where the gravitational potential would then be 
proportional to G^/^. 



VI. DISCUSSION 

In recent papers [l-4|, nonlocality has been introduced into classical gravitation theory 
via a scalar kernel k. However, observational data can provide information about its re- 
ciprocal kernel q. This is similar to the situation in general relativity, where gravitation 
is identified with spacetime curvature, but observations generally do not directly measure 
the curvature of spacetime, except possibly in relativistic gravity gradiometry. We make a 
beginning in this paper in the treatment of the inverse problem of nonlocal gravity. The 
scalar nonlocal kernel k must be determined from observational data that involve the recip- 
rocal kernel q. Our preliminary study involves the Newtonian regime, where the nonlocally 
modified Poisson's equation is investigated in its linearized convolution form. We present a 
detailed mathematical analysis of the resulting Fredholm integral equation using the Fourier 
transform method and prove the existence of the nonlocal convolution kernel A;(u) when 
its reciprocal g(u) satisfies certain physically reasonable conditions. Simple explicit exam- 
ples are worked out in connection with the linear gravitational potential of spiral galaxies. 
To extend our treatment beyond the Newtonian domain, it would be necessary to consider 
relativistic generalizations of the Kuhn kernel along the lines indicated in Sec. Ill of Ref. j^. 

Next, we present a general treatment of the nonlocal and nonlinear modification of Pois- 
son's equation that represents nonlocal gravity in the Newtonian regime. The method of 
successive approximations is then employed to provide a formal solution. The utility of this 
general approach is illustrated for the determination of the gravitational potential of a point 
mass when nonlinearities are assumed to be relatively small. In this case, we recover anew 
the Tohline-Kuhn phenomenological modified gravity approach to the dark matter problem 



in astrophysics 



12 



To place our work in the proper context, we note that nonlocal special relativity (devel- 
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oped since 1993, cf. jsl) and the principle of equivalence imply the necessity of a nonlocal 
generalization of Einstein's theory of gravitation. Here nonlocality is encoded in a nonlocal 



"constitutive" kernel k that must be determined from observation. 



n working out the phys- 



ical consequences of nonlocal gravity, it was soon discovered that it reproduces the 

1980s Tohline-Kuhn phenomenological approach to dark matter as modified gravity 12l. Il3|. 
This connection is the most fundamental contact of the new theory with observation and 
indicates to us that we are on the right physical track. To verify this, we must compute the 
nonlocal kernel k from the rotation curves of spiral galaxies and show that it has the proper 
physical properties expected of such a kernel. Our present paper accomplishes this task. 
That is, we extend the Kuhn kernel q analytically to all space and then use the result to 
solve the inverse problem of finding kernel k by means of Fourier integral transforms. The 
resulting k has indeed just the expected properties and puts the nonlocal theory of gravity 
on a more solid observational foundation. 
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Appendix A: Radial Convolution Kernels 

The purpose of this appendix is to present some useful relations between the radial 
convolution kernels that we employ in Sec. III. We use dimensionless quantities throughout. 
A radial convolution kernel q{u) can be expressed in terms of its Fourier transform g(^) 

as 



?M = ^ / m)M^u)d^- (Al) 



1 

Eqs. f lAip and ( 130|1 form a pair of Fourier sine transforms such that 







1 r°° 

?(0) = 7rT/ (A2) 



2vr2 J, 



CO 



q{0) = Att / r'q{r)dr. (A3) 
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It is interesting to note that Eqs. (133|) and (lA2p imply that 

It is clear from a comparison of Eqs. (1A3P and (149|) that g(0) ~ Md/M. Therefore, we can 
conclude from the discussion in Sec. Ill that g(0) ^ 2/a in Figures 1 and 3, in agreement 
with our numerical results. 

For the two particular examples considered in Sec. Ill, g(0) = (1 + aa)/{47ia'^) for the 
first example given by Eqs. f p^ and fl35|) . and g(0) = 00 for the second example that has a 
central cusp and is given by Eq. ( 15T]) . Let us note that in either case g(^) is square integrable 
over the whole ^ space, so that C,'^q^{C,) is integrable over the radial coordinate : — ?■ 00; 
therefore, the right-hand side of Eq. (lASp is finite. It then follows from Eq. (lASp that k{0) is 
finite in the first example due to the finiteness of q{0), while in the second example, g(0) = 00 
and hence —k{0) = 00, in agreement with the numerical results of Figures 2 and 4. 

Finally, let q'{u) and q"{u) represent respectively the reciprocal kernels given in Sec. Ill 
in the first and second examples; then, 

/// N // N a 1 + a(a + u) ^„ , 

q"{u) - q\u) = , \ . A6 

47r u[a + uY 

Moreover, we find from Eqs. f lA6P and (!30!) that 

s^o \-\ci + T) a + rJ 



It then follows from the lemma given at the end of Sec. II that the right-hand side of Eq. (I A 711 
is positive. Thus for any ^ > 0, q"{Cj > ^'(Oi moreover, as ^ ^ 00, q"{^) — q'{^) — )■ 0, in 
accordance with the Riemann-Lebesgue lemma. 
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